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The force on a macroscopic polarizable body in an inhomogenous electromagnetic field is cal- 
culated for three simple exactly solvable situations. Comparing different approaches we pinpoint 
possible pitfalls and resolve recent confusion about the force density in ferrofluids. 

1 Introduction 

The force on a macroscopic electrically or magnetically polarizable body in an external 
electromagnetic field is one of the most basic and at the same time most relevant problems 
in classical field theory. Although all texts 1-4 on the subject introduce Maxwell stress and 
the force on a (microscopic) dipol there is still considerable confusion on how to use these 
concepts in concrete situations. That these problems are not confined to the classroom is 
exemplified by a recent controversy among experts on the magnetic force density in ferrofluids 

5,6 

In order to discuss the possible problems which may arise when calculating the force on a 
polarizable body it is very useful to study simple and exactly solvable cases. Since a non-zero 
force requires an inhomogenous external field such examples are not abound. The purpose 
of the present paper is to present three simple situations in thermodynamic equilibrium for 
which fields, stresses and forces can all be calculated exactly, using only basic Maxwell theory, 
and to discuss for these situations different methods to calculate the resulting ponderomotive 
forces. The first example will be dealt with in detail whereas for the remaining two only the 
main results will be given. 

2 Basic equations 

Neglecting effects of electro- or magnetostriction the Maxwell stress tensor for a macroscopic 
polarizable medium in equilibrium is of the form 2 ' 4 

T i:i = -i(e E 2 + a*oH 2 )% + E i D j + HiBj, (1) 

where as usual E denotes the electric field, D the electric displacement, H the magnetic 
field, and B the magnetic induction. From this stress tensor one finds the Kelvin expression 

fi : djTij = {PjdrfEi + no(Mjdj)Hi (2) 

for the ponderomotive force density in the medium. Here P := D — eoE denotes the electric 
polarization and M := B//io — H the magnetization of the medium whereas £o and fio are 
the dielectric constant and magnetic permeability respectively of the vacuum. Finally the 
change of the equilibrium free energy of a polarizable body due to a change of the fields is 
given by 3,4 

5F = - I d 3 r (P ■ (SEo + MoM • <5H ), (3) 
Jv 
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where the integral is over the volume of the body, and Eo and Ho denote the external fields 
in the absence of the polarizable material. 

All three expression are consistent with each other and can therefore be used equally 
well to determine the total force on a macroscopic body. In the following examples we will 
always assume linear constitutive relations, i.e. D = £o£ r E and B = /io/ir-H, since we are 
only concerned with effects for which a possible nonlinearity in this relations is not crucial. 

3 Point charge polarizing a half space 

Let us first consider a point charge Q at a = (0, 0, a) a distance a > away from a dielectric 
half space with relative dielectric constant e r > 1 (see fig.[j]). We will neglect possible 
magnetic effects and study only the electric part of the problem. 




Figure 1: Lines of the electric displacement field D for a point charge Q polarizing a 
dielectric half space with e r = 3. The field is symmetric with respect to rotations around 
the z-axis. 



As is well known 2,1 the electrostatic problem can be solved exactly with the help of 
a mirror charge Q' — ~Q(e r — 1)/(£V + 1) at r = —a and an additional auxiliary charge 
Q" = 2Q/(e r + 1) at r = a. Using the coordinate system of fig.[j] the result for the electric 
field is 
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Our aim is to determine the force F with which the point charge attracts the half space. 
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3.1 The physical picture 

The field of the point charge polarizes the dielectric which results in an induced bulk charge 
density p — — VP and a surface charge density a = P n := P • n where P = Eo(s r — 1)E is the 
polarization of the medium and n denotes the normal vector on the surface pointing outward 
of the medium. In our case n — (0, 0, 1). Since the field in the dielectric is a Coulomb field 
of a point charge (cf.(Q)) we have VP = and no bulk charge density is induced. For the 
surface charge we find from (||) setting z = 



a(x,y) = s {e r - 1) 



Q 2 a 
Atteq e r + 1 (x 2 + y 2 + a 2 ) 3 / 2 ' 



(5) 



Summing up the Coulomb forces between Q and the elements of this surface charge we find 
F = (0, 0, F) with 



F 



1 1 



Q 2 . 

47re e r + 1 4a 2 ' 



(6) 



This is an unambigous and physically sound expression for the force in question. It is exactly 
equal to the Coulomb force between Q and the mirror charge Q'. 

3.2 Integrating the Kelvin force 

For the Kelvin force density (@) we find in the present situation 
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Integrating this expression over the halfspace z < the x and y components again average 
to zero and for the z component we get 



F 
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This is clearly different from (^|): the integral over the Kelvin force density does not yield 
the total force on a polarizable body. This is what was observed also in 5 for another 
example. The reason for this discrepancy is however not that the expressions (0) or (0) 
for the Kelvin force density are incorrect but that the total electromagnetic force on the 
polarizable medium involves both bulk and surface contributions whereas (||) only covers the 
bulk part. The expression for the appropriate surface part is given below. 



3.3 Balance of stresses 

By the definition of stress the total force on a body is given by the integral of the stress 
tensor over the complete surface of the body. Paying attention to the fact that the Maxwell 
stress is non-zero also outside the polarizable medium a calculation of the electromagnetic 
force on a body has to integrate the differences in stress slightly above and below the surface 
of the body. In our particular problem the fields decay sufficiently rapidly to zero for r — > oo 
such that only the integral over the x-y-plane contributes. Using (Q) and (0) we find for the 
relevant component of the Maxwell stress tensor 



lim T zz (x,y,z) = 
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on the medium side and 



lim T zz (x,y,z) 
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on the vacuum side. Integrating the difference between these electric stresses one finds for 
the z component of the force 
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This is again different from the correct expression (Q)! The reason for this renewed failure 
is somewhat more subtle. The point is that we have only considered the electric part of the 
stress tensor in the medium whereas the total force must clearly be related to the surface 
integral of the total stress tensor T*° tal . But in equilibrium we must have djTf? = 
(otherwise there would be momentum transport in the medium) and the electric part (j^) 
is hence counterbalanced by non-electric contributions. If our dielectric is a solid, elastic 
stresses would build up, if it is a fluid a non-trivial pressure field would compensate the 
electric stress. In any case the total internal stress in equilibrium must be divergence free 
and hence cannot contribute to the total force on the body. 

The whole stress on the surface is therefore given by the vacuum part ( |To| ) and its surface 
integral indeed yields the correct result @. Incidentally this also explains why the force 
is equal to the Coulomb force between Q and Q': the field for z > is identical for either 
the polarizable half space or the dielectric replaced by the mirror charge Q' and hence the 
integrals over the vacuum stress must be the same. 

Let us finally note that from the definition fa := djT^j of the Kelvin force density it is 
clear that the volume integral of fi must be identical to the surface integral of TL-. Hence the 
elimination of the contribution from the electric stress inside the medium in (O) can simply 
be accomplished by adding . Equivalently ( |ll| ) is nothing but the surface contribution 
that was missed in (Js|) . Using the behaviour of the fields E and D at the interface between 
media with different dielectric constant it is easy to show 4 that 
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Hence the appropriate surface force density complementing the Kelvin body force density is 
given by 
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(13) 



3.4 Change in free energy 

The simplest and most direct way to determine the total force on a polarizable body in 
equilibrium is to use expression (|^) for the variation in the free energy. Specializing to the 
case without magnetic contributions we note first that the change in the electric field Eo in 
the region of the body due to a displacement of the body by an infinitesimal vector 5r is 
given by SFjq = (Sr -V) Eo. On the other hand the corresponding change in free energy is 
related to the force F via SF = — F ■ Sr by the very definition of the total force F. Using the 
condition V x Eo = fulfilled in thermodynamic equilibrium we therefore find 

F = / d 3 r (P • V) E . (14) 
Jv 
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Using ([!]) to determine P as well as 



Q r — a 

Eo r - -f- 15 

47reo k" — or 



and performing the integral we indeed recover the correct expression 



4 Charged line polarizing a parallel cylinder 

This is again an example from electrostatics. A line at distance a from and parallel to the 
y-axis carries a charge Q per length. A dielectric cylinder of radius R, dielectric constant e r , 
and the y-axis as symmetry axis is polarized by the field of the line charge (see fig||) . We 
want to determine the force with which the cylinder is attracted to the line. 




Figure 2: Lines of the electric displacement field D for a charged line polarizing a parallel 
dielectric cylinder with e r — 3 and radius R = a/2, where a denotes the distance between 
the axis of the cylinder and the charged line. The situation is symmetric with respect to 
translations in the y-direction. 

Due to the translational invariance along the y-axis the problem is essentially two- 
dimensional and accordingly we use the notation r = (x, 0, z) in the following. The electro- 
static problem can again be solved by the method of images. To determine the fields outside 
the dielectric one introduces two additional charged lines with densities Q' = —Q(e r — 
l)/(e r + 1) at a distance a' :— R 2 /a from the y-axis and —Q' at the y-axis respectively. 
The field inside the dielectric is determined by replacing the original line by one with charge 
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density Q" = 2Q/(e r + 1). The resulting electric field is given by 
Q 2 r a 
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(16) 



The field inside the dielectric is again divergence free. The induced surface charge can be 
determined from the normal component of the polarization as in the first example and is 
found to be 



a{9) = 



Q e r 



1 



R — a cos 9 



7r e r + 1 R 2 - 2Ra cos 9 



(17) 



where 9 denotes the angle between the radius vector to the surface element in question and 
the z-axis. Using the two-dimensional Coulomb law and summing up the ^-components of 
the forces between the surface elements and the original charged line we find 
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(18) 



As expected this is exactly the sum of the Coulomb forces between Q, Q' and —Q'. 

The most direct way to calculate this force is to use (||) in the form Since P is 

parallel to Eq we have 



(PV)E = PV^o 
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Integrating this expression over the volume of the dielectric cylinder the x-component aver- 
ages to zero whereas the z-component yields back (fig). 



5 Magnetizable cylinder in the field of a straight wire 

Finally we discuss a simple magnetic example for which the electric part will be neglected. 
Consider a straight wire at distance a from the y-axis carrying a current I. A cylinder of 
radius R and magnetic permeability jjL r with the y-axis as symmetry axis is placed into the 
magnetic field of the wire. We want to determine the force per length on the cylinder. 

Due to the well known similarities between the magnetic fields of parallel, stationary 
currents and two-dimensional electrostatics the algebra for this example is very similar to 
the previous one. We introduce two mirror currents, V = I([i r — l)/(/i r + 1) parallel to the 
original one at a distance a' = R 2 /a from the y-axis and —I' along the y-axis. The field 
inside the medium is given by replacing I by I" = 21 /(p r + 1). The magnetic field is then 
found to be 



H(r) = 



I 2 e y x (r — a) 
2ir fi r + 1 \r — a\ 2 
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(20) 
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Figure 3: Lines of the magnetic induction B for a current-carrying straight wire magnetizing 
a parallel cylinder with magnetic permeability fi r — 3. The geometry is identical with fig.^. 

where e y denotes the unit vector in y-direction. 

To clarify the physical picture we note that the magnetization M = (/i r — 1)H correspon- 
des to induced bulk currents of density j b = V x M and induced surface currents j,=Mx n. 
In our case V x M = and 



I H r — 1 a cos 9 — R 
7r fi r + 1 R 2 - 2Ra cos 9 + a 2 



?.W = :^ n D -„ fl , ^ ( 21 ) 



which is analogous to (|T^) . Correspondingly summing up the magnetic forces between these 
surface currents and the original one we find 

F= ^llt^± f (22) 
2vr fi r + 1 a(a 2 - R 2 ) V ' 

This result is most easily derived by using (|J) which in the present case gives rise to 

F = fiQ I d 3 r (M • V) H . (23) 
Jv 

Exploiting the fact that M and Ho are parallel it is easy to reproduce (p2]). 

6 Summary 

Three simple examples have been presented for which the equilibrium force on a polarizable 
body in an external electromagnetic field can be determined exactly. Different ways to do 
this have been elucidated. A first possibility is to integrate the Maxwell stress tensor (Q) 
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over a surface outside the body and completely enclosing it. The second method integrates 
the volume force density (Q) and adds the surface integral over the difference in normal 
stresses (P^/(2£q) + involving the normal components of the polarization and the 

magnetization. The third and in our opinion most simple and straightforward method builds 
on the free energy change (||) and uses the formula 

F= f d 3 r [(P • V) E + Mo (M • V) H ] , (24) 
Jv 

where Eo and Ho denote the fields in the absence of the polarizable medium. Note that 
in the case where the fields do not vary much over the volume of the body this expression 
reduces to the well known relation 

F = (P-V)Eo + /io(M-V)H , (25) 

with "P and A4 denoting the total polarization and magnetization respectively of the body. 
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